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Abstract. The 3:2 spin-orbit resonance between the rotational and orbital motions of Mercury (the periods are 
P$ = 56.646 and Pa = 87.969 days respectively) results from a functional dependance of the tidal friction adding 
to a non-zero eccentricity and a permanent asymmetry in the equatorial plane of the planet. The upcoming space 
missions, MESSENGER and BepiColombo with onboard instrumentation capable of measuring the rotational 
parameters stimulate the objective to reach an accurate theory of the rotational motion of Mercury. For obtaining 
the real motion of Mercury, we have used our BJV model of solar system integration including the coupled 
spin-orbit motion of the Moon. This model, expanded in a relativistic framework, had been previously built in 
accordance with the requirements of the Lunar Laser Ranging observational accuracy. We have extended the 
BJV model by generalizing the spin-orbit couplings to the terrestrial planets (Mercury, Venus, Earth, and Mars). 
The updated model is called SONYR (acronym of Spin-Orbit N-BodY Relativistic model). As a consequence, 
the SONYR model gives an accurate simultaneous integration of the spin-orbit motion of Mercury. It permits to 
analyze the different families of rotational librations and identify their causes such as the planetary interactions or 
the parameters involved in the dynamical figure of the planet. The spin-orbit motion of Mercury is characterized 
by two proper frequencies (namely $ = 15.847 years and ^ = 1066 years) and its 3:2 resonance presents a second 
synchronism which can be understood as a spin-orbit secular resonance (II = 278 898 years). A new determination 
of the mean obliquity is proposed in the paper. By using the SONYR model, we find a mean obliquity of 1.6 
arcminutes. This value is consistent with the Cassini state of Mercury. Besides, we identify in the Hermean 
librations the impact of the uncertainty of the greatest principal moment of inertia (C /MR 2 ) on the obliquity 
and on the libration in longitude (2.3 milliarcseconds and 0.45 arcseconds respectively for an increase of 1% on 
the C /MR 2 value). These determinations prove to be suitable for providing constraints on the internal structure 
of Mercury. 
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1. Introduction 

Before 1965, the rotational motion of Mercury was as- 
sumed to be synchronous with its orbital motion. In 
1965, Pettengill & Dyce discovered a 3:2 spin-orbit reso- 
nance state by using Earth-based radar observations (the 
Mercury's rotation period is = 58.646 days while the 
orbital one is P\ = 87.969 days). This surprising reso- 
nance results from a non-zero eccentricity and a perma- 
nent asymmetry in the equatorial plane of the planet. In 
addition, the 3:2 resonance strongly depends on the func- 
tional dependance of the tidal torque on the rate of the 
libration in longitude. Moreover the 3:2 resonance state is 
preserved by the tidal torque (Colombo & Shapiro 1966). 
The main dynamical features of Mercury have been es- 
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tablished during the 1960's in some pioneer works such as 
Colombo (1965), Colombo & Shapiro (1966), Goldreich & 
Peale (1966) and Peale (1969). Goldreich & Peale (1966) 
notably studied the probability of resonance capture and 
showed that the 3:2 ratio is the only possible one for a 
significant probability of capture. In addition, in a tidally 
evolved system, the spin pole is expected to be trapped 
in a Cassini state (Colombo 1966; Peale 1969, 1973). The 
orbital and rotational parameters are indeed matched in 
such a way that the spin pole, the orbit pole, and the so- 
lar system invariable pole remain coplanar while the spin 
and orbital poles precess. The reader may find in Balogh & 
Campieri (2002) a review report on the present knowledge 
of Mercury whose the interest is nowadays renewed by two 
upcoming missions: MESSENGER (NASA, Solomon et al. 
2001) and BepiColombo (ESA, ISAS, Anselin & Scoon 
2001). 
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Our work deals with the physical and dynamical causes 
that contribute to induce librations around an equilib- 
rium state defined by a Cassini state. In order to wholly 
analyze the spin-orbit motion of Mercury and its rota- 
tional librations, we used a gravitational model of the so- 
lar system including the Moon's spin-orbit motion. The 
framework of the model has been previously constructed 
by Bois, Journet & Vokrouhlicky (BJV model) in accor- 
dance with the requirements of Lunar Laser Ranging (LLR 
thereafter) observational accuracy (see for instance a re- 
view report by Bois 2000). The approach of the model 
consists in integrating the A^-body problem on the ba- 
sis of the gravitation description given by the Einstein's 
general relativity theory according to a formalism derived 
from the first post-Newtonian approximation level. The 
model is solved by modular numerical integration and 
controlled in function of the different physical contribu- 
tions and parameters taken into account. We have ex- 
tended this model to the integration of the rotational mo- 
tions of the terrestrial planets (Mercury, Venus, Earth, 
and Mars) including their spin-orbit couplings. The up- 
dated model is then called SONYR (acronym of Spin- 
Orbit 7V-bodY Relativistic model). As a consequence, us- 
ing SONYR, the Af-body problem for the solar system and 
the spin motion of Mercury are simultaneously integrated. 
Consequently we may analyze and identify the different 
families of Hermean rotational librations with the choice 
of the contributions at our disposal. 

Starting with the basic spin-orbit problem according to 
Goldreich & Peale (1966), we have computed a surface of 
section for the Mercury's rotation showing its very regular 
behavior. We have calculated again the proper frequency 
for the spin-orbit resonance state of Mercury. Using our 
model, an important part of the present study deals with 
the main perturbations acting on the spin-orbit motion 
of Mercury such as the gravitational figure of the planet 
as well as the planetary effects and their hierarchy. A de- 
tailed analysis of the resulting rotational librations due to 
these effects is presented and described in the paper. A 
new determination of the Hermean mean obliquity is also 
proposed. Moreover, we identify in the Hermean librations 
the impact of the variation of the greatest principal mo- 
ment of inertia on the instantaneous obliquity and on the 
libration in longitude. Such a signature gives noticeable 
constraints on the internal structure of Mercury. 

2. Geometry of the spin-orbit coupling problem 

According to Goldreich & Peale (1966), we consider the 
spin-orbit motion of Mercury with its spin axis normal to 
the orbital plane. The orbit is assumed to be fixed and un- 
variable (semi-major axis a and its eccentricity e). The po- 
sition of Mercury is determined by its instantaneous radius 
r while its rotational orientation is specified by the angle 
9. The orbital longitude is specified by the true anomaly 
/ while the angle 9 — f measures the angle between the 
axis of least moment of inertia of Mercury and the Sun- 
to- Mercury line (see Fig.rj). According to these assump- 



tions, the dynamical problem of the spin-orbit motion of 
Mercury is reduced to a one-dimensional pendulum-like 
equation as follows : 

o CM 

Ce + ^(B-A)^^ S m2(6-f) = (1) 

where G is the gravitational constant, M Q the solar mass, 
and A, B, and C the principal moments of inertia of 
Mercury. This equation of motion has only a single degree 
of freedom for the spin-orbit coupling, the characteristic 
angle of rotation 9, but depends explicitly on time through 
the distance r to the planet and the non-uniform Keplerian 
motion of the true anomaly /. As a consequence, it is a 
problem not reducible by quadrature and the equation 
is non-integrable (Wisdom 1987). 

^ Planet A e "/ 
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Sun 


/ Peri apse 



Fig. 1. Geometry of the spin-orbit coupling problem, r is 
the radius vector, / the true anomaly and 9 is the char- 
acteristic angle of rotation. 



In order to know the structure of the phase space of 
the Mercury's rotation, a surface of section for its spin- 
orbit coupling is very useful. Let be a — \J ^ B q the 
asphericity of the Mercury's dynamical figure combining 
the principal moments of inertia A, B, and C. The equa- 
tion |TjQ) becomes : 

d+\n 2 (^) 3 a 2 sin 2(0- /) = (2) 

where n is the orbital mean motion, / and 9 as de- 
fined above. The numerical integration of the equation J3J 
yields to the Figure showing stroboscopically, i.e. once 
point per orbit at the periapse, various trajectories of the 
phase portrait of Mercury (9 and l/nd9/dt on the axes 
x and y respectively). Figure J21a) illustrates the different 
quasiperiodic librations around the 1:2, 1:1 (synchronous), 
3:2, and 2:1 rotation states while Figure J5|b) shows a 
zoom of the present 3:2 resonance state for Mercury. Let 
us notice that in the 1:2 rotation state, the two libration 
areas are shifted by 7r/2 with respect to the other rotation 
states. The former rotation rate of Mercury was higher 
and the planet could have experienced large and chaotic 
variations in obliquity at some time of the past (Laskar & 
Robutel 1993). Due to its slowing down by tidal despin, 
the Mercury's rotation has then encountered different res- 
onance states, from higher orders to the one trapped in 
the present state. Goldreich & Peale (1966) have shown 
that this present state (the 3:2 resonance) was the first 
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Fig. 2. Surface of section of the Mercury's spin-orbit coupling (a — 0.0187, e = 0.206). The rate of change of the 
orientation is plotted versus the orientation of Mercury defined by the angle 9 at every pericenter passage. The spin 
axis is constrained to be normal to the orbital plane. The plot (a) makes in evidence the different quasiperiodic 
librations around the 1:2, 1:1 (synchronous), 3:2, and 2:1 rotation states. The plot (b) is a zoom of the real 3:2 
resonance state of Mercury. The chaotic zone is microscopic and non visible at this scale. 



one with a substantial probability of capture. This one 
indeed depends crucially on the functional dependence of 
the tidal torque acting on the spin-orbit motion. This de- 
pendence is expressed by the variation rate of the longi- 
tude libration angle 7 as defined below. In the end, the 
stabilization of Mercury at the 3:2 spin-orbit resonance is 
due to permanent asymmetry on the equatorial plane, as 
well as its non-zero eccentricity equal to 0.206 (Colombo 
1966, Colombo & Shapiro 1966). 

According to the Chirikov resonance overlap criterion 
(1979), the chaotic behavior appears when the asphericity 
of the body is larger than the following critical value : 



1 



(3) 



where e is the orbital eccentricity. In the case of Mercury, 
a — 0.0187 is lower than the critical value a cr — 0.2701. 
As a consequence, the spin-orbit behavior is regular. The 
zoom in Figure b) shows indeed that the separatrix 
surrounding the 3:2 resonance state is very small. The 
width of the associated chaotic zone is then estimated to 
the order of 10 -43 (Wisdom et al. 1984). Consequently, 
tidal friction pulls Mercury accross the chaotic separatrix 
in a single libration period. 

From the equation ijljl. it is possible to obtain an in- 
tegrable approximated equation using the spin-orbit res- 
onance, the spin rate 6 being commensurable with the 
mean orbital motion n. Following Murray & Dermott 



(2000), by introducing a new variable 7 



pn where 



p parametrizes the resonance ratio (p = | in the case 
of Mercury), one may expand the equation in form- like 
Poisson series. Taking into account that 7 n, one av- 
erages all the terms over one orbital period, and finally 
obtain the following equation : 



3 2 
7 +-n 



B-A 
C 



H{p, e) sin 2 7 = 



(4) 



where H {p, e) is a power series in eccentricity. In the case 
of Mercury, this expression is written as follows : 

7 123 , 



(5) 



2 16 

Finally, the proper frequency of the Mercury's spin-orbit 
motion is : 



UJ 



,B-A 



H(p,e)\ 



(6) 



which by using the values listed in Table|21gives the proper 
period of 15.830 years. 

Balogh & Giamperi (2002) developed the equation J2J 
and obtained the following expression : 



7 + a G 2 o 9 (e) sin[2 7 + (1 - q)M] = 



(7) 



where the G2o q coefficients are eccentricity functions de- 
fined by Kaula (1966). The G201 coefficient is equal to 
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ff(f,e). a = = 1.76 10~ 4 for the Mercury case. 

Figure |3| presents then the numerical solution of the for- 
mula (thanks to G. Giampieri, private communica- 
tion). The angle 7 describes a periodic behavior with a 
period equal to the revolution of Mercury. The amplitude 
of 42 arcseconds (as) depends on the value of oto. The be- 
havior of the 7 angle is nicely matched by the approximate 
formula of Jhen & Corral (2003) : 

7 = 0o sin M + fa sin 2M (8) 

where fa = a (G 2 oo - G 202 ) and fa = ^f(G 20 -i - G 203 ). 
The authors noticed that the ratio of the two amplitudes 
{(j) and fa) does not depend on the a parameter, and 
K = = -9.483. 




50 100 150 200 250 



t (days) 

Fig. 3. The spin-orbit solution of Mercury in the planar 
case (equation[7} plotted over 250 days. Arcseconds are on 
the vertical axis and days on the horizontal axis. Short- 
term librations have a period of 87.969 days (the orbital 
period of Mercury) and 42 as of amplitude. 



However, the above equations describe a spin-orbit mo- 
tion of Mercury where the spin axis is normal to the orbital 
plane while the orbital motion is Keplerian. 

3. The extended BJV model : SONYR 

3.1. Theory 

In order to wholly analyze the spin-orbit motion of 
Mercury and its rotational librations, we have enlarged 
a gravitational model (called BJV) of the solar system in- 
cluding the Moon's spin-orbit motion. The accurate the- 
ory of the Moon's spin-orbit motion, related to this BJV 
model, was constructed by Bois, Journet & Vokrouhlicky 
in accordance with the high accuracy of the LLR obser- 
vations (see previous papers : Bois et al. 1992; Bois & 
Journet 1993; Bois & Vokrouhlicky 1995; Bois et al. 1996; 
Bois & Girard 1999). The approach of the BJV model con- 
sists in integrating the A^-body problem (including trans- 
lational and rotational motions) on the basis of the grav- 



itation description given by the Einstein's general rela- 
tivity theory. The equations have been developped in the 
DSX formalism presented in a series of papers by Damour, 
Soffel & Xu (Damour et al. 1991, 1992, 1993, 1994). For 
purposes of celestial mechanics, to our knowledge, it is 
the most suitable formulation of the post-Newtonian (PN) 
theory of motion for a system of N arbitrary extended, 
weakly self-graviting, rotating and deformable bodies in 
mutual interactions. The DSX formalism, derived from the 
first post-Newtonian approximation level, gives the post- 
Newtonian representation of the translational motions of 
the bodies as well as their rotational ones with respect to 
the locally transported frames with the bodies. 

Gravitational fields of the extended bodies are parame- 
terized in multipole moment expansions : (M£, S^) define 
the mass and spin Blanchet-Damour multipoles charac- 
terizing the PN gravitational field of the extended bodies 
while (G^,Hl) are tidal gravitoelectric and gravitomag- 
netic PN fields. Because we do not dispose of dynami- 
cal equations for the quadrupole moments M^ b , and al- 
though the notion of rigidity faces conceptual problems 
in the theory of relativity, we have adopted the 'rigid- 
multipole' model of the extended bodies as known from 
the Newtonian approach. Practically this is acceptable 
since the relativistic quadrupole contributions are very 
small. Consequently and because it is conventional in geo- 
dynamical research to use spherical harmonics analysis 
of the gravitational fields with the corresponding notion 
of harmonic coefficients (G^S^J, the quadrupole mo- 
ments MA have been expressed in those terms, according 
to reasons and assumptions given in Bois & Vokrouhlicky 
(1995). Gravitational figures as well as the figure-figure 
interactions of the bodies are then represented by expan- 
sions in spherical harmonics (Borderies 1978, Shutz 1981). 
Moreover, internal structures of solid deformable bodies, 
homogeneous or with core-mantle interfaces, are repre- 
sented by several terms and parameters arising from tidal 
deformations of the bodies (both elastic and anelastic). 
More details and references on these topics are given in 
the above quoted papers related to our works concerning 
the theory of the Moon's spin-orbit motion. 

The BJV model, as described above, has been ex- 
tended to the spin-orbit integration of the terrestrial 
planets (Mercury, Venus, Earth, and Mars). This new 
model is henceforth called SONYR (for Spin-Orbit N- 
bodY Relativistic model). In the present paper framework, 
the SONYR model is devoted to the detailed analysis of 
the complete spin-orbit motion of Mercury. 

The simultaneous integration of the solar system, in- 
cluding the Mercury's spin-orbit motion, uses a global 
reference system given by the solar system barycenter. 
Nevertheless, let us recall that local dynamically non- 
rotating frames show a slow (de Sitter) rotation with re- 
spect to the kinematically non-rotating frames. As a con- 
sequence, the reference frame for the Mercury's rotation 
is affected by a slow precession of its axes transported 
with the translational motion of Mercury. In the Earth's 
case, the de Sitter secular precession of the Earth reference 



N. Rambaux and E. Bois: Theory of the Mercury's spin-orbit motion and analysis of its main librations 



5 



frame is close to 1.92 as/cy (see Fukushima 1991, Bizouard 
et al. 1992, Bois & Vokrouhlicky 1995). Consequently, the 
real rotation of Mercury has not to be expressed in an 
incrtial system fixed in space, but in a local dynamically 
non-rotating frame fallen down in the gravitational field of 
the Sun. Because of the proximity of Mercury to the Sun, 
its de Sitter precession may be expected quite significant. 

In the end, the SONYR model and its analysis method 
take into account (i) the experience in post-Newtonian 
gravitation in the definition of reference frames required to 
deal with rotational motions combined with translational 
ones, and (ii) the modern knowledge of dynamical systems 
for studiing librations as quasi-periodic solutions accord- 
ing to the axiomatic presented in Bois (1995). We can state 
that the model is not Newtonian but rather 'Newtonian- 
like', resulting from truncation of the fully post-Newtonian 
(DSX) framework. In the present paper, we deal with the 
Newtonian- like librations (classical physical librations), 
while the formally relativistic contributions (relativistic 
librations and de Sitter precession of the Mercury's refer- 
ence frame) will be analyzed in a forthcoming paper. 

3.2. Method 

The model is solved by modular numerical integration and 
controlled in function of the different physical contribu- 
tions and parameters taken into account. The iV-body 
problem (for the translational motions), the rotational 
motions, the figure-figure and tidal interactions between 
the required bodies are simultaneously integrated with the 
choice of the contributions and truncations at our disposal. 
For instance, the upper limits of the extended figure ex- 
pansions and mutual interactions may be chosen as fol- 
lows : up to I = 5 in the Moon case, 4 for the Earth, 2 for 
the Sun while only the Earth-Moon quadrupole-octupole 
interaction is taken into account (see previous papers). 
The model has been especially built to favor a systematic 
analysis of all the effects and contributions. In particular, 
it permits the separation of various families of librations 
in the rotational motions of the bodies. 

The non-linearity features of the differential equations, 
the degree of correlation of the studied effect with respect 
to its neighbors (in the Fourier space) and the spin-orbit 
resonances (in the Moon and Mercury's cases), make it 
hardly possible to speak about 'pure' effects with their 
proper behavior (even after fitting of the initial condi- 
tions). The effects are not absolutely de-correlated but 
relatively isolated. However, the used technique (modu- 
lar and controlled numerical integration, differentiation 
method, mean least-squares and frequency analysis) gives 
the right qualitative behavior of an effect and a good quan- 
tification for this effect relative to its neighbors. In the 
case of the particular status of the purely relativistic ef- 
fetcs, their quantitative behaviors are beyond the scope of 
the present paper and will be discussed in a forthcoming 
work. When a rotational effect is simply periodic, a fit of 
the initial conditions for a set of given parameters only 



refines without really changing the effect's behavior. The 
amplitudes of librations plotted in Figures ^] and ^| are 
then slightly upper bounds. 

The precision of the model is related to the one re- 
quired by the theory of the Moon. One of the aims in build- 
ing the BJV model (at present included inside SONYR) 
was to take into account all phenomena up to the precision 
level resulting from the LLR data (i.e. at least 1 cm for 
the Earth-Moon distance, 1 milliarcsecond (mas) for the 
librations). For reasons of consistency, several phenomena 
capable of producing effects of at least 0.1 mas had been 
also modeled (the resulting libration may be at the ob- 
servational accuracy level). Moreover, in order to justify 
consistence of the Moon's theory, this one had been ad- 
justed to the JPL ephemeris on the first 1.5 years up to a 
level of a few centimeter residuals. In the other hand, the 
internal precision of the model is only limited by the nu- 
merical accuracy of the integration. Thus, in order to avoid 
numerical divergence at the level of our tests for Mercury, 
computations have been performed in quadrupole preci- 
sion (32 significant figures, integration at a 10 -14 internal 
tolerance) . 

3.3. Terminology 

In order to de-correlate the different librations of Mercury, 
we use the terminology proposed in Bois (1995), which 
is suitable for a general and comparative classification of 
rotational motions of the celestial solid bodies. This ter- 
minology derives from a necessary re-arrangement of the 
lunar libration families due to both progress in the Moon's 
motion observations (LLR) and modern knowledge of dy- 
namical systems. 

Traditionally, the libration mode called physical libra- 
tion is split up according to the conventional dualism 
"forced- free". The forced physical librations are generally 
related to gravitational causes while the free librations 
would be departures of the angular position from an equi- 
librium state. These cuttings out contain ambiguities and 
redundancies discussed in previous papers (Bois 1995, Bois 
2000). Formally, the free librations are periodic solutions 
of a dynamical system artificially integrable (by a conven- 
tion of writing related to specific rates of the spin-orbit 
resonance, for instance 1:1), whereas the forced librations 
express, in space phase, quasi-periodic solutions around a 
fixed point (the system is no longer integrable). Moreover, 
any stable perturbed rotation of celestial solid body con- 
tains imbricate librations of different nature, and those are 
too strongly overlapped to keep the traditional classifica- 
tion. 

In the present terminology, the libration nature, its 
cause and its designation are linked up. Two great libra- 
tion families serve to define the physical librations, namely 
the potential librations and the kinetical librations. They 
simply correspond to a variation energy, potential or ki- 
netical respectively. For libration sub-classes, the designa- 
tion method is extensive to any identified mechanism (see 
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Table 1. Our initial conditions at 07.01.1969 (equinox J2000). (a) Mean values derived from the SONYR model; (b) 
Seidelmann et al. (2002) 





Mercury 






Rotation angles 




ip = 48.386 deg 


ip = 0.0 deg/day < ip > = 


-0.616 10~ 7 deg/day (a) 


9 = 7.005 deg 


6 = 0.0 deg/day < 6 > = 


-0.267 10~ 8 deg/day (a) 


tp = 299.070 deg 


ip = 6.138505 deg/day (b) 770 = 


1.6 amin (a) 



Table 2. Parameters of Mercury, (a) JPL; (b) Anderson et al. (1987); (c) Milani et al. 2001. 



Mercury 
Bulk quantities 



Mass (GMq) = 4.9125. 10~ 10 (a) 

Equatorial radius (km) = 2439 (b) 

J 2 = (6.0±2.0).10~ 5 (b) 

C22 = (1.0±0.5).HT 5 (6) 

C/MR 2 = 0.34 (c) 



Table 3. Our results for the spin-orbit motion of Mercury. The spin-orbit period verifies the relation : P = 2P\ = ZPu, 



Mercury 
Spin-orbit characteristic period 



$ (I s * proper frequency) — 


15.847 years 


(2 nd proper frequency) = 


1066 years 


II (orbital precession) = 


278 898 years 


Pa (orbital period) = 


87.969 days 


P v (rotational period) = 


58.646 days 


P (spin-orbit period) = 


176.1 days 



more details in Bois 1995). The terminology permits eas- 
ily the separation of various families (see the Moon's case 
described in a set of previous papers). These librations are 
called direct when they are produced by torques acting on 
the body's rotation. They are called indirect when they 
are produced by perturbations acting on the orbital mo- 
tion of the body. Indirect librations derive from spin-orbit 
couplings. 

A specificity of the SONYR model with its method 
of analysis is to isolate the signature of a given pertur- 
bation. The SONYR model allows indeed the identifica- 
tion of relationships between causes and effects including 
interactions between physics and dynamics, such as the 
dynamical signature of a core-mantle interaction (called 
centrifugal librations) . 

3.4. Parameters and initial conditions 

In the computations presented in the paper, the required 
dynamical parameters and general initial conditions come 
from the JPL DE405 ephemeris (Standish 1998). However, 
concerning the parameters related to the Mercury's ro- 
tation (second-degree spherical harmonics C20 and C22), 
which are not included in the JPL ephemeris, our model 



uses those given by Anderson et al. (1987) (see Table 0). 
Besides, up to now it does not exist any ephemeris of 
the Mercury's rotation. As a consequence, to build ini- 
tial conditions for the Hermean rotation (described by an 
Eulerian sequence of angles ip, 9 and tp defined below in 
Section 4-1), we use the following principle : Assuming 
the polar axis of Mercury normal to its orbital plane, 
we obtain ip = Q and 9 = i where f2 and i are respec- 
tively the ascending node and the inclination of the orbit 
of Mercury on the Earth equatorial plane (which is the 
reference frame used in the DE405 ephemeris). The long 
axis of Mercury being pointed towards the Sun at its peri- 
apse allows to fix the ip angle of polar rotation. The value 
of tp is found in Seidelmann et al. (2002). We use at last 
ip = and 9 — 0; these two variables reach to mean values 
generated by the complete spin-orbit problem of Mercury : 
— 0.616 x 10~ 7 deg/day and — 0.267 x 10~ 8 deg/day respec- 
tively. The numerical integrations presented in the paper 
start from these initial conditions related to the planar 
problem for Mercury; they are listed in Tabled Departure 
from the planar case is understood as the integration of 
physics included in SONYR. 

In the other hand, for the computations carried out 
in this paper, the global reference frame O'X'Y'Z' is 
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Fig. 4. The rotational motion of Mercury expressed in the ecliptic reference frame OXYZ by using the 3-1-3 Eulerian 
sequence (if), 9, ip), and plotted over 3000 years. Degrees are on the vertical axis and years on the horizontal axis. The 
value of 278 898 years for the secular part of the if) angle is related to the planetary interactions. The period of 1066 
years in the three angles is the libration of the spin-orbit secular variable. The period of 15.847 years in the ip angle 
is the proper frequency in longitude of Mercury. 



given by a reference system centered on the solar system 
barycenter, fixed on the ecliptic plane, and oriented at 
the equinox J2000. The rotational motion of Mercury is 
evaluated from a coordinate axis system centered on the 
Mercury's center of mass Oxyz relative to a local dynami- 
cally non-rotating reference frame, OXYZ, whose axes are 
initially co-linear to those of O' X'Y' Z' . In the framework 
of the present paper without purely relativistic contribu- 
tions, let us note that axes of OXYZ remain parallel to 
those of O'X'Y'Z'. 

The TV-body problem for the planets of the solar sys- 
tem and the Mercury's spin-orbit motion are simultane- 
ously integrated in the SONYR model. Concerning the 
rotational equations written in a relativistic framework, 
the reader may refer to Bois & Vokrouhlicky (1995). In a 
Newtonian approach, these equations amount to the clas- 
sical Euler-Liouville equations of the solid rotation (see 
e.g. Goldstein 1981). We follow the formalism and the 
axiomatic expanded in Bois & Journet (1993) and Bois 
(1995) for the definition of the different rotational ele- 
ments as well as the used terminology. Let us simply pre- 
cise that 1 is the angular momentum expressed in Oxyz 
and is related to the instantaneous rotation vector u) as 
follows : 

1=00" (9) 
where (I) is the tensor of inertia for the body. According 
to classical assumptions, (/) is reduced to three princi- 
pal moment of inertia A, B, and C . The gravity field of 



Mercury is essentially unknown. The tracking data from 
the three fly-by of Mariner 10 in 1974-75 have been re- 
analyzed by Anderson et al. (1987) to give a low accuracy 
estimate of the normalized coefficients C20 and C22 (see 
Table values are expressed in the body-fixed frame of 
the principal axes of inertia). The principal moments of 
inertia A and B are then infered from C20 and C22 by the 
following formulae (Ferrari et al. 1980) : 

1 C 20 - 2C 22 + -JL, (10) 



MR 2 



B 



a 



MR 2 ~ 20 ' ' MR 2 

Let us note that the parameter C /MR 2 is not constrained 
within these relations. We present in Section the vari- 
ations oiC/MR 2 by analyzing its signature into the libra- 
tions of Mercury. 

4. The librations of Mercury 

4.1. Planetary perturbations 

For Mercury assumed to be a rigid body reduced to three 
oblateness coefficients, the general expression for a torque 
coming from a disturbing point mass m at the vectorial 
position ru (r is the instantaneous distance) is written as 
follow : 



2C : 



MR 2 
C 



22 



N 2 = 



3GV 



- U X (I) u 

rJ 



(11) 
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Fig. 5. Mercury's spin-orbit secular resonance. The top 
panel presents the evolution of ip an d ^ both including 
a secular term. The bottom panel presents the behavior 
of £ = ip — plotted over 3000 years. Let us underline 
that this £ angle does not present any secular term, show- 
ing then a synchronism between both angles, ip and £7. 
On the top pannel, degrees are on the vertical axis while 
arcminutes are on the bottom panel; years are on both 
horizontal axes. 



The ip angle expresses the nodal precession of the equa- 
torial plane of Mercury with respect to the ecliptic plane. 
It splits up in a periodic term with a period "J = 1066 
years and a secular one II = 278 898 years. \& is the second 
proper frequency of the Mercury's spin-orbit coupling. It 
can be analytically approximated by the following formula 
(used in the Earth's case by Goldstein 1981) : 



. ,'3n 2 C- A 

27:/ h- z — cost 



(12) 



This analytical expression is suitable for an axis symmet- 
ric body. It is not the case of Mercury. However, assuming 
Mercury as a symmetric top rotating about its smallest 
axis of inertia, with an average equatorial moment of iner- 
tia A = A + B , one finds for this period 1300 years (by using 
the values given in Tabled- The difference between the 
two values permits to appreciate the departure of Mercury 
from a symmetric body. In the other hand, the dynamical 
behavior of ip coming from SONYR is due to the direct 
effect of the solar torque by the way of the true dynamical 
figure of Mercury. In the 2-body problem, Sun- Mercury, 
^ = 1066 years and II = years. Related to the planetary 
interactions, II = 278 898 years expresses the departure 
from the 2-body problem. 



where a is the mean distance between the two bodies. 
Such a torque due to the Sun is the dominant one acting 
on the rotational motion of Mercury. The solar torque is 
indeed responsible of the global dynamical behavior of the 
rotation of Mercury. Due to Venus, the resulting torque 
is about of the order of 10 5 times smaller, that is quite 
negligible. 

Figure presents the rotational motion of Mercury 
including only the solar torque in the rotational equa- 
tions and taking into account simultaneously the whole 
A^-body problem for the Sun and the planets (the plane- 
tary interactions inducing indirect effects on the rotation 
of Mercury). In this Figure, the Euler angles ip, 6, ip related 
to the 3-1-3 angular sequence describe the evolution of the 
body-fixed axes Oxyz with respect to the axes of the local 
reference frame OXYZ. Let us recall the definition used 
for these angles : ip is the precession angle of the polar 
axis Oz around the reference axis OZ , 9 is the nutation 
angle representing the inclination of Oz with respect to 
OZ, and p is the rotation around Oz and conventionally 
understood as the rotation of the greatest energy (it is 
generally called the proper rotation). The axis of inertia 
around which is applied the proper rotation is called the 
axis of figure and defines the North pole of the rotation 
(Bois 1992). Let us remark that in the Figure0](in other 
Figures involving ip as well) plotted over 3000 years, we 
have removed the mean rotation of 58.646 days in the (f 
angle in order to better distinguish the librations. We may 
then clearly identify the first proper frequency of 15.847 
years (to be compared to the analytical determination, 
namely 15.830 years given in Section 2). 
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Fig. 6. The rotational motion of Mercury expressed in 
the ecliptic reference frame OXYZ by using the 3-1-3 
Eulerian sequence (ip,9,ip), and plotted over 500 days. 
Degrees are on the vertical axes and days on the horizontal 
axis. The periods of 58.646 and 87.969 days in the ip and 
ip angles respectively, express precisely the 3:2 spin-orbit 
resonance rate. The period P=176.1 days appears clearly 
in these plots. 



The dynamical behaviors of ip and Q (the ascending 
node of the orbit) are quite superimposed as shown in 
Figure03(top panel). As it is mentioned by Beletski (1986), 
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Fig. 7. Interactions of the planets on the Hermean rotational motion by the way of the spin-orbit couplings over 3000 
years. Degrees are on the vertical axes and years on the horizontal axis for the three panels. In the black line case, 
the problem is reduced to the Sun and Mercury. In the broken line case, the interactions between the Sun, Mercury 
and Venus are taken into account. With the dots line the later case includes Jupiter in addition. The whole planetary 
interactions are integrated in the cross line case (except for Pluto) . The dots line case defines our called standard case 
sufficient for preserving the 3:2 resonance ratio. 



a second synchronism is generally involved in a Cassini 
state. However, in the Mercury's case, it is not a periodic 
synchronism as it is the case for the Moon. Using SONYR, 
we make easily in evidence the periodic 18.6 year synchro- 
nism in the lunar spin-orbit resonance while it is about a 
secular behavior in the Hermean spin-orbit resonance. The 
3:2 resonance of Mercury introduces a mechanism of spin- 
orbit secular resonance qualitatively analogous to the or- 
bital secular resonances combined with mean motion reso- 
nances (as for instance in the 2:1 case, see Bois et al. 2003). 
As it is shown in Figure [S] (bottom panel), the spin-orbit 
secular resonance variable £ = ip — ft does not present 
any secular term. As a consequence, ip and O on average 
precess at the same rate equal to n, confirming then the 
mechanism of spin-orbit secular resonance. We find that 
£ librates with the particular frequency of 1066 years. 

In addition, in order to give a detailed inspection of the 
short periods involved in the rotational motion of Mercury, 
Figure El presents the solution plotted over 500 days. The 
P v rotation period of 58.646 days appears in the ip an d 
9 angles. Whereas the mean rotation of 58.646 days is 
removed in the ip angle (as in Figure 0J) , the signature 
of the P\ orbital period of 87.969 days is clearly visible 
(this angle is called libration in longitude of 88 days in 
literature). A third period appears in the ip and 9 angles, 
namely 176.1 days. This one results from the 3:2 spin-orbit 
resonance (P = 2P\ = 3P V ). 



Figured presents the planetary interactions acting on 
the rotational motion of Mercury by the way of its spin- 
orbit couplings (i.e. indirect effects of the planets on the 
Mercury's rotation). In the black line case the problem 
is reduced to the Sun and Mercury. In this 2-body prob- 
lem, the orbital plane does not precess as it is clear in the 
ip and 9 angles without secular terms. The secular vari- 
ations rise up from the departure of the 2-body problem 
(as it is visible with the broken, dots, and cross line cases 
in Figure 01. In the broken line case, the interactions be- 
tween the Sun, Mercury and Venus are taken into account. 
With the dots line, the later case includes Jupiter in ad- 
dition. The whole planetary interactions are integrated in 
the cross line case (except for Pluto). We show that Venus 
is the planet which induces the greatest secular term. After 
Venus, the role of Jupiter is dominant, and this 4-body 
problem (Sun, Mercury, Venus, and Jupiter) defines our 
"standard" case used in our Section 4-3 for the analysis 
of the Hermean librations. The rate of secular variations 
in the Mercury's rotation between all planetary interac- 
tions (cross lines) and our standard case (dots lines) is 
11.8 as/cy (as : arcseconds) in the 9 nutation angle and 
1.9 amin/cy (amin : arcminutes) in the ip precession angle. 
These values should be used as corrective terms in analyt- 
ical theories of the rotational motion of Mercury. Let us 
emphasize that the spin-orbit motion of Mercury coming 
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from our standard case is sufficient for preserving the 3:2 
resonance ratio between the two modes of motion. 

Starting with the initial conditions defined in Section 
3.4 (where in particular the initial obliquity of Mercury is 
equal to zero), the SONYR model permits obtaining the 
dynamical behavior of the Hermean obliquity by its si- 
multaneous spin-orbit integration. The variables i, 9, f2, ip 
from SONYR substituted in the following relation : 



cos 77 



cos 1 cos f 



sin 1 sm 



' cos(fi — ip) 



(13) 



produce the instantaneous obliquity 77 plotted over 3000 
years as presented in Figurc[S] Such a behavior for 77 gives 
easily a mean obliquity 770 = 1.6 amin. We show at present 
that this mean value is quite consistent with the Cassini 
state of Mercury. Let be the following equation established 
by Peale (1988) and reformulated by Wu et al. (1997) and 
coming from the Cassini laws : 



MR 2 
C 



sin (i + 770) 



n sin?7o[(l + cos77o)G 2 oiC2 2 - cos 770G210C20] 



(14) 



where fj, is the precessional angular velocity of the 
Hermean orbit, namely 27r/II years^ 1 , while 



G 



7 



2111 



123 
Iff* 



and G 



210 



are eccentricity functions defined by Kaula (1966). i is the 
inclination of the orbital plane of Mercury relative to a 
reference system precessing with the orbit. This inclina- 
tion varies between from 5° to 10° while the eccentricity 
varies from 0.11 to 0.24, over 10 6 years (Peale 1988). As a 
consequence, 770 obtained with belongs to [1.33, 2.65] 
amin. Let us note that a conventional value of 7 amin is 
often given in literature. Such a value, outside the inter- 
val of possible values, is very probably incorrect as already 
claimed by Wu et al. (1997). 

4.2. Principal figure librations 

4.2.1. Signature of the C/MR 2 coefficient on the 
rotational motion of Mercury 

Let us consider at present the disturbing torques acting on 
the rotational motion of Mercury and as a consequence in- 
ducing direct librations. This section focuses on the libra- 
tions related to the dynamical figure of the planet. Such 
librations are called principal figure librations (Bois 1995). 
We assume the Sun reduced to a point mass while the 
gravity field of Mercury is expanded in spherical harmon- 
ics up to the degree 2. We express the solar torque acting 
on the figure of Mercury according to the equation (|11|) . 

The first coefficients of the Hermean gravity field have 
been determined with the Mariner 10 probe (Anderson 
et al. 1987). We use these values for C 2 o — —Ji and C22 
given in Tabled In order to complete the Hermean tensor 
of inertia (coefficients A, B, C), the C/MR 2 principal mo- 
ment of inertia is required (see equation Ijlljl l. Its value is 
related to the internal density distribution of the planet 
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Fig. 8. Dynamical behavior of the r\ instantaneous obliq- 
uity of Mercury plotted over 3000 years. Arcminutes are 
on the vertical axis and years on the horizontal axis. Over 
this time of integration, the behavior is simply described 
by a period of 1066 years and 3.2 amin of amplitude. The 
770 mean obliquity is 1.6 amin. 



according to the polar axis of Mercury (rotation of great- 
est energy about the smallest principal axis of inertia). 
For an homogeneous planet, such a normalized value is 
equal to 0.4. We use a nominal value of 0.34 (Table (2J 
used by Milani et al. (2001) and coming from an internal 
structure model of Mercury including three layers (crust, 
mantle and core). 
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Fig. 9. Rotational behavior of Mercury for two different 
values of its greatest principal moment of inertia. Degrees 
are on the vertical axes and years on the horizontal axes. 
The computation with C/MR 2 =0A is plotted through the 
dashed lines; this value corresponds to an homogeneous 
planet. The broad lines are obtained with C/MR 2 =Q.3A; 
this value derives from an internal structure model of 
Mercury including three layers. On the bottom panel (</? 
angle), the dashed lines are shifted away 1 degree in order 
to distinguish them from the broad lines. 
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Fig. 11. Signature of the C /MR 2 coefficient on the obliquity of Mercury for two different values shifted from 1%. 
Arcminutes are on the vertical axis of the top panel while milliarcseconds are on the vertical axis of the bottom panel; 
days are on the horizontal axis for both panels. The maximal amplitude of the signature of C /MR 2 on the obliquity is 
2.39 mas. On the top panel, the dashed lines are shifted away 0.001 amin in order to distinguish them from the broad 
lines. 
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Fig. 10. The signature of the libration of Mercury when 
the principal moment of inertia along the axis of figure 
varies from 1 % in the top panel (a), to 2% in the mean 
panel and 15 % in the bottom panel (c). Degrees are on 
the vertical axes and years on the horizontal axes. 



Figure presents the rotational behavior of Mercury 
computed over 10 000 years in our standard case with 
two different values of its greatest principal moment of 
inertia : (i) the computation with C/MR 2 =QA is plot- 
ted with the dashed lines (homogeneous planet); (ii) the 
broad lines are obtained with C/Mi? 2 =0.34 (three layer 



model of Mercury). On the bottom panel (ip angle), the 
dashed lines are shifted away 1 degree in order to distin- 
guish them from the broad lines. This Figure[!|] (especially 
the angle) shows how the C/MR 2 coefficient value and 
the constant of precession 'J? (2 nd proper frequency) are 
related. With C/MR 2 =0A, W = 1254.01 years while with 
C/Afi? 2 =0.34, W = 1066 years. Besides, C/MR 2 and $ 
(I s * proper frequency) are also linked and for evaluating 
this relation, Figure ITU1 shows the variations Aip obtained 
by differentiation : (i) on the top panel AC = 1%, (ii) 
on the middle panel AC = 2%, (iii) on the bottom panel 
AC = 15% (0.34+ 15%(0.34) = 0.4). The beats are sig- 
natures related to the variations in 

4.2.2. Signature of the C/MR 2 coefficient on the 

obliquity and on the libration angle in longitude 

One of the main objectives of the BepiColombo and 
MESSENGER missions is to measure the rotation state of 
Mercury, up to an accuracy allowing to constrain the size 
and physical state of the planet's core (Milani et al. 2001, 
Solomon et al. 2001). Consequently, the two missions have 
to determine the four following parameters : C20, C22, ?7, 
and ip that are sufficient to determine the size and state of 
the Mercury's core (see Peale 1988, 1997). Combining C20, 
C22, and 77, one obtains the C/MR 2 coefficient while with 
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Fig. 12. Signature of the C '/MR 2 coefficient on the Hermean libration angle in longitude for two different values 
shifted from 1%. Degrees are on the vertical axis of the top panel and arcseconds are on the vertical axis of the bottom 
panel; days are on the horizontal axis for both panels. The maximal amplitude of the signature of C /MR 2 on the 
libration angle in longitude is 0.45 as. 
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C-22 and tp, one obtains C m /MR 2 (i.e. the C/MR 2 coef- 
ficient for the mantle) . The validity condition of the first 
combination is that the dynamical behavior of the core 
has to follow the one of the mantle over a period of time 
at least the one of II (assertion 1). The validity condition 
of the second combination is that the dynamical behavior 
of the core has not to be coupled to the one of the man- 
tle over a period of 88 days (P\) (assertion 2). These two 
conditions linked together imply some constraints on the 
nature of the core-mantle interface (Peale 1997). In order 
to reach such an objective, the BepiColombo mission has 
to obtain a value on the C/MR 2 coefficient with an ac- 
curacy of 0.003, i.e. 1% and therefore foresees measuring 
the libration angle and the obliquity with an accuracy of 
3.2 and 3.7 as respectively (Milani et al. 2001). 

Our SONYR model gives (i) the true relation between 
the three parameters (C/MR 2 , r\, ip), and (ii) the up- 
per bounds of the impact of C m /MR 2 on the ip angle. 
FigureirTlpresents the impact of the C/MR 2 coefficient on 
the instantaneous obliquity rj. In these plots (Figures 1111 
and I12|) , the spin-orbit motion of Mercury is integrated 
within the whole solar system with an initial obliquity of 
1.6 amin, which is the mean obliquity of Mercury evalu- 
ated in Section J^.l. The top panel of Figure ITU expresses 
the dynamical evolution of r\ computed over 500 days with 
C/MR 2 =0.M (black lines) and C/M i? 2 =0.3434 (dashed 
lines). Dashed lines are shifted from 0.001 amin in order 



to distinguish the two different kinds of lines. Figure ITTI 
shows also how the instantaneous obliquity of Mercury 
differs from its 1.6 amin nominal value. The bottom panel 
shows by differentiation the signature of the 1% variation 
of C/MR 2 on 77. The maximal amplitude crest to crest is 
of the order of 2.3 mas within the characteristic period of 
176.1 days related to the 3:2 ratio of the Mercury's spin- 
orbit resonance. 

Figure H*2l presents the signature of the C/MR 2 coeffi- 
cient on the tp libration angle in longitude. The top panel 
expresses the behavior of 99 computed over 500 days with 
C/Mi? 2 =0.34 (black lines) and C/M i? 2 =0.3434 (dashed 
lines). One may compare this FigurelT^lto the Figurc[3]re- 
sulting from the usual analytical resolution of the Eulerian 
equation 0] (thanks to G. Giampieri, private communica- 
tion). Let us note that the angle 7 defined in Figure is 
equal to the angle <p plotted in Figure H*2l The later only 
gives a simple double sine curve with an amplitude of 42 as 
while the solution of the SONYR model includes the cou- 
plings between the three rotational variables as well as the 
indirect couplings due to planetary interactions (we notice 
that in the two FigureslT2landl3*lthe amplitude of libration 
is of the order of 40 as). Let us note that the FigurclT2lcor- 
responds to the libration related to the C/MR 2 coefficient 
of the planet without core-mantle couplings. Let us add 
that in Peale (1972), the amplitude of tp is related to the 
C m coefficient by assuming that the assertion 2 quoted up- 
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per is true. On the contrary, our first results on this topic 
make in evidence the existence of a faint coupling. This 
core-mantle coupling will be presented in a forthcoming 
paper. The bottom panel of Figure IT21 shows by differen- 
tiation the signature of the 1% variation of C/MR 2 on 
(p. The maximal amplitude within the period of about 88 
days (i.e. the signature of Pa) is of the order of 0.45 as. 

In conclusion, signatures of the indeterminacy of 1% in 
C/MR 2 on the obliquity and on the libration in longitude 
are 2.4 mas and 0.45 as respectively. What is very faint 
(may be too much) with respect to the expected accuracy 
forecasted in the BepiColombo mission. 

4.3. variing the obliquity 

Because the initial obliquity value is unknown, we test in 
this last section the impact of the indeterminacy of this 
value on the spin-orbit motion of Mercury. The results 
are presented in Figure El plotted over 3000 years; top 
panel: the effect on the nutation angle 9, middle panel: 
the effect on the orbital inclination i, and bottom panel: 
the effect on the instantaneous obliquity r\. On each panel, 
three curves are related to three different initial values of 
rj, namely amin (black lines), 1 amin (dashed lines), and 
2 amin (dot lines) . In the bottom panel, the amplitudes of 
these librations are of the order of 1.4 amin with a period 
of 1066 years. For any initial value of 77 S [0, 3.2] amin, 
the mean value of rj, let be 770, is equal to 1.6 amin, which 
is in good agreement with the determination of 7/0 in a 
previous section. We may claim that 770 = 1 ■ 6 amin. 

For obtaining such a mean obliquity by measurements, 
let us underline that the theoretical behavior of 77 points 
out to fit the observations by a sine function taking into 
account the long period = 1066 years with an amplitude 
of 1.6 amin. 

5. conclusion 

The 3:2 spin-orbit resonance between the rotational and 
orbital motions of Mercury results from a functional de- 
pendance of the tidal friction adding to a non-zero ec- 
centricity and a permanent asymmetry in the equato- 
rial plane of the planet. The upcoming space missions, 
MESSENGER and BepiColombo with onboard instru- 
mentation capable of measuring the rotational parameters 
stimulate the objective to reach an accurate theory of the 
rotational motion of Mercury. 

Starting from our BJV relativistic model of solar sys- 
tem integration including the coupled spin-orbit motion of 
the Moon, we have obtained a model generalizing the spin- 
orbit couplings to the terrestrial planets (Mercury, Venus, 
Earth, and Mars). The updated model is called SONYR 
(acronym of Spin-Orbit N-BodY Relativistic model). It 
permits to analyze and identify the different families of 
rotational librations. This work has been carried out for 
Mercury in the present paper. 

The spin-orbit motion of Mercury is characterized by 
two proper frequencies (namely $ = 15.847 and \P = 1066 
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Fig. 13. Impact of the initial obliquity on the nutation 
angle 9 (top panel), on the orbital inclination i (mid- 
dle panel), and on the instantaneous obliquity rj (bottom 
panel). Degrees are on the vertical axes of the top and 
middle panels while arcminutes are on the vertical axis of 
the bottom panel; years are on all horizontal axes. On each 
panel, three curves plotted over 3000 years correspond to 
three different initial values of 77 : amin (black lines), 1 
amin (dashed lines), and 2 amin (dot lines). 



years) and its 3:2 resonance presents a second synchro- 
nism which can be understood as a spin-orbit secular res- 
onance, (n = 278 898 years). A new determination of the 
mean obliquity has been proposed in the paper. By using 
the SONYR model, we have found a mean obliquity of 
1.6 amin. This value is consistent with the Cassini state 
of Mercury. Besides, we have identified in the Hermean li- 
brations the impact of the uncertainty of the greatest prin- 
cipal moment of inertia (C/MR 2 ) on the obliquity and on 
the libration in longitude (2.3 mas and 0.45 as respectively 
for an increase of 1% on the C/MR 2 value). These deter- 
minations prove to be suitable for providing constraints on 
the internal structure of Mercury. The direct core-mantle 
interactions will be presented in a forthcoming paper. 
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Appendix A: free rotation 

Let us assumed Mercury isolated in space; in this sense, 
its rotation is free and the Euler-Liouville equations for its 
rotation are written without right hand side, i.e. without 
any external disturbing torques. If we add the assumption 
of a rigid Mercury, we are in the Euler-Poinsot motion 
case (whose solutions are the well-known Eulerian oscilla- 
tions). Without explicitly integrating such equations, the 
assurance of integrability in the Poincare sense can be ob- 
tained by some theoretical simple considerations. Indeed, 
whatever being the triplet of generalized coordinates used 
to describe the spatial attitude of a solid body in a fixed 
frame, one knows that there exists four independant inte- 
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grals of motion : the hamiltonian H , and the three compo- 
nents Lx , Ly , Lz of the angular momentum (in OXYZ) . 
Four integrals of motion for three degrees of freedom, the 
problem is then integrable and even over-integrable. One 
does not lose the generality of the problem choosing for 
instance A < B < C . The choice A < B < C makes pos- 
sible to write the general solution of the system under a 
form involving the elliptical functions of Jacobi (Landau 
& Lifchitz 1969). By convention, let us adopt that the 
resulting oscillations in space be called the Eulerian oscil- 
lations, expressing exclusively the oscillations of the non- 
perturbed rotation of the rigid body. From this resolution, 
we obtain the Eulerian frequencies : 



and 



lL 2 ~2AH 
C(C - A) 



= n, 



(A.i) 



(A.2) 



that give for Mercury the two periods fi/ and a (where 
L is the angular momentum in OXYZ, H the energy), a 
and (3 are the dynamical coefficients of the body related 
to A, B and C by the following relations : 



and [3 ■ 



.4 



B 



(A.3) 



With the previous initial conditions given in Table ^ 
(ip = 8 = 0), the free rotation of Mercury is reduced to 
an elementary rotation of the tp angle only (J7/ = lu z ). 
Starting with the physical mean values of tp, 8, and 770 
evaluated in the present paper, flf = 58.646 days and 
a = 964.88 years. The period of 58.646 days corresponds 
to the Hermean polar rotation while the one of 964.88 
years means the global period of rotation of Mercury in 
R 3 , which is equivalent to the Earth's classical Euler pe- 
riod, namely 305 days (Lambeck 1980). In the case of the 
Moon, this period is equal to 148.129 years. 

Using the SONYR model reduced to the free rotation 
of Mercury, we obtain the components of the instanta- 
neous rotation vector u> in the body- fixed system Oxyz, 
as presented in Figure lA~Tl u z is well found constant while 
the Euler period is equal to 964.92 years. 
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